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Abstract 



We demonstrate that the nuclear colhsion geometry (i.e. impact parameter) can be determined 
with 1.5 fm accuracy in an event-by-event analysis by measuring the transverse energy flow in the 
pseudorapidity region 3 < I77I < 5 with a minimal dependence on collision dynamics details at the 
LHC energy scale. Using the HIJING model we have illustrated our calculation by a simulation of 
events of nucleus-nucleus interactions at the c.m.s energy from 1 up to 5.5 TeV per nucleon and 
various type of nuclei. 



The measurement of the colhsion impact parameter is a very important practical problem of 
relativistic heavy ion physics. A part of future experimental programmes on the LHC accelerator 
will be devoted to relativistic heavy ion collisions Significant efforts will be focused on the 
establishment of the fundamental laws and quite general rules of nucleus-nucleus interactions, the 
discovery of a new state of QCD matter, i.e. quark-gluon plasma (QGP), and the study of the 
properties of strongly exited nuclear matter p. It is expected that QGP might be produced in the 
central nucleus-nucleus collisions at extremely high energy density eo ~ 1 GeV/ fiv? (at the LHC 
energy scale eo ~ 1 TeV / frn?). 

On the other hand, for studies of diffractive phenomena, properties of a coherent pomeron and 
collective nuclear effects § it is necessary to select peripheral events with a large collision impact 
parameter. Assuming that the collision impact parameter is measured the experimentally observed 
effects can be compared with theoreticall predictions of expected signals of a "new" physics: parton 
energy losses in nuclear matter, monojet to dijet ratio, quarkonia suppression, correlated jet and 
W^,Z^ production and so on. 

In the present paper we demonstrate a method of estimation of a collision impact parameter in 
the event-by-event analysis. 

As the collision impact parameter can not be measured in experiments directly, it would be 
necessary to find an experimentally measurable 6-dependent variable. We suggest to use the total 
transverse energy Et (total energy E) produced in the one event of the nucleus-nucleus collision for 
the collision impact parameter estimation. 

We will show that there is a correlation between the transverse energy and the impact param- 
eter, and will present our numerical calculations performed on the basis of the HIJING model 
at the RHIC and LHC energy scales. We argue that the measurement of Et produced in the for- 
ward direction (large (pseudo)rapidity values) allows one to avoid possible uncertainties in the b 
determination. 

The multiplicity and transverse energy production in the nucleus-nucleus collisions in the ultra- 
relativistic energy domain is considered as a combination of hard processes with p-r > Po and soft 
particle production. A transverse energy flow produced in hard processes described by perturbative 
QCD is associated in the main with minijet production, i.e. jets with ~ 4 GeV 0. 

The average transverse energy carried by (mini)jets in the rapidity intervals Ay is related to the 
collision impact parameter b by the formula: 

{Erib, y/SNN,P0,Ay)) = TAA{b)ajet{y/^NN,Po)Ay{ET^), (1) 

where (E^^) is the average transverse energy per a (mini)jet in Ay interval, ajet is the cross-section 
of minijet production in the parton model at the pp- level. The differential distribution da/dEx is: 
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{y/sNN,Po,^y) = 9^ / dp'^dyidy2 Xifi{xi,p'^)x2fj{x2,PT) 
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dt dt 1 + Oki 

where po is the pQCD cut-off parameter, xi and X2 are the fractional momenta of the initial par- 
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tons i and j, and yi and y2 are the rapidities of outgoing partons, da'^^^'^^ / di is the parton-parton 
scattering cross section. The summation runs over all parton species and the factor K k, 2 is used 
to correct the lowest order pQCD rates for the effects of next-to-leading order terms. The value of 
crjet{y/sNN,Po)Ay{E^^) then is: 

(Tjet{y/SNN,Po)Ay{Ep) = dpldyidy2 ^ Xifi{Xi,p^)x2fj{x2,PT) X 

[—^(5, t, u) + — ^(S, t)]-^- (3) 
dt dt 1 + Oki 

Experimentally known effects of modification of quark and gluon structure functions by nuclear 
medium called parton shadowing @] has not been included here. To take into account the shadowing 
effect we should modify the formula by multiplying parton distributions by a corresponding 
correction term: 

fiixi,PT)fjix2,PT) ' Rtixi,pl)fiixi,PT)Rfix2,PT)fjix2,pl), (4) 

here the ratio Rf^j = fi,j/Aix,p'^)/ fij/N{x,p'^), where fi,j/N{x^p%) is the parton structure function 
for a free nucleon and fi,j/A{x,p^) is the corresponding parton distribution in a proton inside the 
nucleus. 

The nuclear density overlap function of two colliding nuclei at a given impact parameter T^Ai^) 
calculated in the assumption of the Wood-Saxon nuclear density distribution pAir): 

TAA{\h\) = j d'fTA{f^TA{b-r}, (5) 
where r is a 2-dimensional vector defining the interaction point. The nuclear thickness function is 



TA{\r\) = j dzpAi\J\f\' + z^). (6) 

The expression (|l|) relating and b includes a term arisen from hard processes with pt >2 GeV 
only. For a more correct estimation of the collision impact parameter we should take into account in 
addition the part of a total transverse energy flow produced in soft interactions. 



{Err ft ^ TAA{b)asofUnci. {E^'Y°f\ (7) 

where {E^'^)'^"^* is the averaged transverse energy per particle produced by soft interactions. But soft 
processes can not be calculated by pQCD applications and we should to use some phenomenological 
models 0, |^ for the estimation of a soft part of the total energy flow. 

We think that a permissible pseudorapidity region for an the measuring energy is limited by 
the following reason. Main signatures of a possible QGP production are expected in the central 
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(pseudo)rapidity region. One of the discussed features of such a state of nuclear matter is energy 
losses of scattered partons in final state interactions with a dense nuclear matter called jet quenching 
0. Among other effects originated by jet quenching one may expect a significant modification in the 
distributions of the transverse energy flow and charged multiplicity dET/di], dE^/dr], and drich/dr] 
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Indeed, an indication is found for ultrarelativistic energy domain concerning the appearance of a 
wide bump in the interval —2<r]<2 over a pseudorapidity plateau of such distributions due to jet 
quenching. Fig.l demonstrates the evolution of the effect with collision impact parameter variation. 
It is interesting to note that a secondary interaction effect such as a jet quenching modifies only the 
central rapidity part, while both pseudorapidity regions (3 < |?7| < 5) remain practically unchanged. 
Therefore the large rj region can be used for the collision impact parameter estimation with minimal 
dependence on possible signals of new physics in the central (pseudo)rapidity region. 

The total cross-section of AA collisions has been calculated in the framework of a HIJING hybrid 
model of nucleus-nucleus interactions [Q, where the cross-section of hard processes has been defined 
by the formula (^). The contribution of the soft part of a produced particle spectrum, i.e. small 
PT-processes has been simulated by using the FRITIOF and DPM models [0, §]. In these models 
hadrons are considered as relativistic strings exited at hadron interactions. Calculations shows that 
for central PbPb collisions at the LHC energy the inclusive cross-section of soft interactions at the 
pp-level is equal to 57.0 mb while the inclusive cross-section of hard processes is equal to 54.3 mb. 
The parton shadowing effect has been taken into account. 

However, it is obvious that the average transverse energy of partons produced in hard processes 
{Et)-''^'^ (3 -i- 5 GeV for \ri\ < 0.5) is larger than the average transverse energy of soft partons 
{EtY°^* ~ 0.4 GeV. This fact reduces more strongly the relative contribution of soft processes in the 
total transverse energy production. In the LHC energy domain already 80 % of the total transverse 
energy is calculated by perturbative QCD application. This allows one to reduce ambiguities of Et 
calculations induced by the use of phenomenology models to take into account small p^-processes. 

In the framework of the HIJING model we have simulated 10.000 events of PbPb, NbNb, CaCa 
interactions at 5.5 TeV per nucleon. The dependence of the total transverse energy produced in the 
pseudorapidity interval 3 < |?7| < 5 on the collision impact parameter is presented in fig.0. The 
variable Et is connected with b. Some fluctuations in nucleus-nucleus collision dynamics limit the 
impact parameter estimation to a 1.5 fm precision. 

The correlation of the same type is presented in fig.|^ with variation of a collision energy. It is 
shown that the increase of the collision energy leads to an improvement of the accuracy the impact 
parameter definition. 

The correlation curve for the total energy flow is of the some shape. The bulk of the energy 
produced in the forward direction run up to 10 ^ 100 TeV . This allows one to measure the total 
energy with high accuracy and to reduce experimental errors for the b estimation. 

We should remark that uncertainties are associated with the use of various parton shadowing 
models and the sets of structure functions lead to an ambiguity in Et — b correlation. In more detail 
the influence of these aspects on the total (transverse) energy flow is discussed in ref [TT . 
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Figure 1: Differential distribution of the total transverse energy dEx/dri (GeV) over pseudorapidity 77 
for 10000 minimum bias PbPb collisions at -s/smn— 5.5 TeV/nucleon with various impact parameter. 
Normalized per number of events. 
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Figure 2: Correlation between the transverse energy flow per collision Et (GeV) in the pseudorapidity 
direction (3< \ri\ <5) and collision impact parameter b (fm). From top to bottom: PbPb, NbNb, 
CaCa collisions at ^/sNN=5.5 TeV/nucleon. 
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Figure 3: Correlation between the transverse energy flow per collision Et (GeV) in the pseudorapidity 
direction (3< |?7| <5) and collision impact parameter b (fm) for PbPb collision at -y/sNN—d.b, 3, 1, 
0.5 TeV/nucleon. 
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